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DISCLAIMER 



Although each program has been tested by its contributor, 
no warranty, express or implied, is made by the contributor 
*or any User's Group, as to the accuracy and functioning of 
the program and related program material, nor shall the fact 
of distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or any User's Group, in connection 
therewith. 



COMMON USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter, ink or pencil) 

Program No. Date 

Program Name: 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

Comment 



2. Does the program do what the abstract says? Yes No^ 

Comment 

3. Is the description clear, understandable, and adequate? Yes Nc> 

Comment 



4. Are the Operating Instructions understandable and in sufficient detail? Yes No^ 

Comment _ 

Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable ? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable)? Yes No_ 

Comment 



6. Does the object program run satisfactorily? Yes No_ 

Comment 



7. Number of test cases ru n Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No> 

Comment 



8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment 

9. Were all necessary parts of the program received? Yes No_ 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 

Company 

Address 



IBM Corporation Company 
Program Information Department 
40 Saw Mill River Road 
Hawthorne, New York 10532 _ 
Attn: PREP FORM COORDINATOR US6rS Gr0up C0Cle 



THIS REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW AND 
EVALUATION PROCEDURE. NONME MBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVALUATION. 
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POLTJOXIAL RC0RES8I0K PROGRAM 
FOB 1HB I. B. M. CARD l620 



0. DTJCBTRA, JR. 
General Foods Research Center 
555 Sooth Broadway 
Tarrytown, K. T. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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DE3CRIPTI0H 

This prograa obtains by least squares the coefficients of a polyuria! 
y m b^I^b^x-x) 1 , vhere 1 may range frost 1 to at most 9, for as many as 
250 pairs (y, x). Prograa switches provide the option of obtaining predic- 
tions y at each of the data points, and the option of obtaining predictions 
for a specified range of x's. Output may be either on the console typewriter 
(fixed point) or from the card punch (floating point). Either or both x and 
y may be transformed. Logarithm, exponential, and reciprocal transformations 
are included. 

CARP PREPARATION 

Data Input consists of 1 header card, H cards each containing a y and an 
x, and, if predictions are called for, a card containing the lowest x for 
which predictions are wanted, an increment, and the highest x. 

HEADER CARD (NO FLAGS) 

Cols. 1-3 
Col. k - 5 
Col. 6 



Col. 7-8 
Col. 9-10 



Col. 11 - 12 
Col. 13 - 1U 

DATA CARDS (MO FLAGS) 

Either x or y may appear first on the card. Column 6 of the HSADKR CARD 
will specify the order. Both must be entered in fixed point form. The format 
is described in the attached description of the CARD C0R7SRSI01 SUBMUUTUR. 



xxx - number of pairs ( x, y) » H 
xx •* highest order wanted - 01 through 09 
x » x, y reversal code 

if y precedes x 

1 if x precedes y 

xx ■ transformation code for y 
xx ■ transformation code for x 

00 * no transformation 

01 ■ replace y or x with log 10 jr or log^ Q x 

02 » replace y or x with e y or e X 

03 - replace y or x with l/y or l/x 
xx ■ order at which residuals should start 
xx - order at which predictions should start 
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VmXCTXM MTOItlOi CARD (WO FLAGS) 

The format requirement for the prediction deflnltloa card it the saw 
as for the data cards. Tb» sequence on the card la (l) the lowest x for which 
predictions are wanted, (2) the lncreaent for x, and (5) the highest x for 
vtaich predictions are wanted. 

These x's will he In the sane units as the original x inputs. All x's 
will be transformed In the sane manner. 

SWITCH SETT DIGS 

With all four switches in the off position, the output consists only of 
the coefficients, their associated suns of squares, and the F ratios using 
the residual scan square as the estimate of error. These F ratios change as 
the residual mean square changes. This output will he in fixed point form 
(6 decimals) from the console typewriter. This output will be punched in 
floating point form If Switch ) is on. 

If Switch 1 Is on and if the order of the polynomial Is not less than 
the digits In columns 11 and 12 of the HEADER CARD, the residuals program 
will be entered. The output will consist of the run identification (i.e., 
the sequence of entry), the transformed y, the corresponding prediction y, 
and the residual, y - y, di Tided by the standard deviation of an observation 
at the x Input. The output will be from the typewriter in fixed point form 
if Switch k Is off and from the card punch In floating point form if Switch k 
is on (in tills case, the transformed x is also Included). 

If Switch 2 Is on and If the order of the polynomial is not less than the 
digits In columns 1} and 1% of the HEADER CARD, the output will Include predic- 
tions for x's specified by the Prediction Definition Card. This output medium 
Is controlled by Switch 4, as for the residuals. The output will consist of 
the untrans formed x, the prediction, and the standard deviation of an observa- 
tion at the specified x. 

All switch settings may be changed during the running of the program. 
If Switch 2 Is off Immediately following the data entry, the Prediction De- 
finition Card will not be read. If at some subsequent time predictions are 
wanted, the instruction "READ PREDICTIOB DWDimORS ROW" will be typed and 
the card will then be read. 



e 
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COHPUTATIOMAL PROCEDURE 

It is necessary to find the coefficients of the polynomial of order D 
given by y ■ b o +EbjX*, where 1 ranges from 1 to D. A disadvantage in such 
a solution is that the computational accuracy becomes increasingly worse, 
because of the loss of significant leading digits. In this program an 
alternative but equivalent polynomial, y - b^Eb^x-x) 1 Is fit to the data, 
In order to decrease the loss in accuracy. 

The equations to be solved are of the form 

Ey " V -^Kx-x) ♦ V**-*) 2 4 • • • + V<*- 5 ) D 

Ey(x-x) - bE(x-x) ♦b 1 E(x-x) 2 +b 2 E(x-x) 5 ♦ . . . +b J) E(x-x) D * 1 

5Hx.i) 2 -b o l(x-x) 2 +b 1 E(x-x) 3 +b 2 E(x-x) ,f + . . . ♦b D E(x-x) D+2 



Ey<x-x) D -h Kx-x) D ♦b 1 E(x-x) D+1 +b 2 E(x-x) D+2 + . . . ♦b D E(x-x) 2D 

In order to determine the error terms it is also necessary to have E(y-y) . 
All summations above extend from 1 to >. E(x-x) is, of course, zero, within 
rounding error. 

The data are transformed (if transformations are specified) and then 
stored in memory. At the same time Ex and Ey are being obtained. A state- 
ment of the analysis being performed is printed. If the logarithm of y and 
the reciprocal of x are specified, the statement "LOGARITHM OF Y VERSUS 
RECIPROCAL OF X PLUS" will appear, followed by the average of the reciprocals 
Of x. 

The sums of the powers of x-x are obtained up to the 20 power, and the 
sums of the products of y and powers of x-x are obtained up to the D power. 
E(y-y) 2 is also obtained. 

The matrix of sums of powers is obtained in triangular form (the matrix 
is symmetric), and a 1 x 1 inverse is formed. The inverse is increased in slse 
as the order of the polynomial is Increased. The inverse matrix la obtained 
for use in the residuals and predictions parts of the program. 

The inversion procedure follows the bordering method described by 
Faddeeva (Computational Methods of Linear Algebra, by T. I. Faddeeva, Dover 
Publications). The method is described below, with some notations! changes. 



© 
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The system of equations can be written in matrix form aa X*t « X *XB. 
Letting A - X*X, C - X% and D - A" 1 , the solution Is B - A"V» DC * 

Let the subscript n denote the site of the matrices at some point In 
the computation. Then 



n-1 

D n-1 
r n-l 



n-1 



'n-1 



B 



n-1 



'n-1 



In this notation A and D are matrices of size n by n. 
n n 

A , and D , are matrices of size n-1 by n-1, U , and r . are vectors of 
n-i n-i ' n-x n-x 

site n-1, B and C are vectors of size n. B . 

' n n n-i. 

n- 1 . •»* A nn' d nn' V 

The t—i i j required by this method consists of the matrix A Q in 

triangular form, the vector B Q , and an auxiliary vector Y Q ^ of size 

n-1, and a fev temporary locations. The equations to form the new inverse 



are vectors of size n, B , and C , are vectors of size 
' n-i n-i 

and C are scalars. 
n 



' nn n 



nJ. ~*n-lfr 



a nn " T n-1 u n-l 

- -d ▼ , 

nn n-1 



"n-1 



-1 
A n-1 



n-1 -n-1 



The first equation forma the auxiliary vector, the second forms the nev 
diagonal, the third provides a nev rov (or column) for the expanded ln- 



Tfae vector 



so that D then is the A for 
n 



verse, and the fourth equation corrects the previous Inverse, 
the incr eas e d n. 

The equations for the coefficients are not given by Faadeeva, but are 
found to bet 
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b - d (C - B' u .) 
n nn * n n-1 n-1 7 

^n " b n/ 4 nn " ^n-l^n-lVl' 

MV B n-1 " P»Tious B n _ 1 ♦ (C Q - B^u^r^ 

The first equation gives the new coefficient, the second gives the sum of 
squares attributable to this coefficient, and the third equation corrects 
the previous coefficients. 

RESIDUALS AND PREDICTIONS 

The coefficients provide a basis for obtaining predictions for specified 
x's, whether they be the input x's or those defined by the Prediction De- 
finition Card. The equation is y m *> +*?l m £> i (x o -x) i ' , vhere x q is the value of 
x at which the prediction is made. 

The prediction equation can be written in matrix form as y»XjB, where 

B is the vector of coefficients, B - [b b....b-.] and X^ is the vector of 

' p o 1 D o 

powers of x -x, X'-[l, x -x. ..(x -x) ]. The estimated variance of the pre- 

o ' o . ' 2° °p 
diction is Xj(X*X) X q s , where s is the residual mean square, which Is the 

only estimate of error at this point. The observations at x Q are assumed to 

be distributed with estimated variance s 2 around the prediction, so that the 

variance of an observation is [l + X'tX^J'^X ]s 2 . The square root of this 

o o 

quantity is the "ERROR OH Y" in the predictions part of this program and the 
divisor for the residual, y-y, in the residuals part. 

OPgRATDIG msTKUCTIONS 

Put the l620 in the manual mode and press Reset . Load the polynomial 
regression program. The Instructions as to the switch settings will be typed, 
and then the rest of the program is loaded. 

Make the switch settings, and then put the header card followed by the 
(y,x) data cards into the reader hopper. Press Reader Start on the 1622 and 
Start on the 1620. 

If Switch 2 is on initially, the Prediction Definition Card must follow 
the last data. If Switch 2 is off initially and turned on after the data 
entry, the computer will halt for the Prediction Definition Card, unless this 
card has been loaded into buffer. 
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After the data entry the average * vill be printed (if transformations 
of the x's ere specified, tale will be the average of the trmnsforaed x'e). 
The 1620 will sees to here stopped after tola, but the iuu of povers end 
croee products are being obtained (the Multiplier neone will be nor log). 

If no transformations are specified, and if aore then 250 pairs (y,x) 
are read, than the program anat be reloaded If transforaatlons are vented In 
the neat polynomial regreeelon. 



© 




EXAMPLE 

An example consisting of 9 pairs of data was set up to illustrate 
the program. The Instructions are illustrated at the top of the following 
page 1. The header, data, and prediction definition cards are listed at the 
bottom of the page. 

The header specifies 9 pairs of data, that at most a cubic is wanted, 
i.e., y = b Q +b 1 (x-x) + b 2 (x-xj^+ bj(x-x)', that the y precedes x on the Input 
cards, that neither y nor x should be transformed, that residuals should 
start with the quadratic, and that predictions should start with the cubic. 
The predictions will be obtained for x » 5, 10, and 15. 

The complete typewriter output Is shown on page 2 and the complete 
card output on page 3. 
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ORDER 


COEFFICIENT 


SUM OF SQUARES 
1*964623000 
c.2f 3224261 


00 


5126888887 


01 


5020895083 
COEFFICIENT 


ORDER 


SUM OF SQUARES 


00 


5126839220 


"7«Tt • J J J 


01 




C.9l37?4?AC. 


02 


471475C124 


4720746232 


RUN 


TRANSFORMED Y 


PREDICTION 


001 


It 1 IOOOOOO 

J 1 p 1 www 


5110481222 


002 


itt QOOAAnn 


ci fL&f flQ3R 


003 


Cf 17finftAflfl 
3 1 1 / uuuuw 




004 




Cf 7QR97ACA 


005 


? f 2400O000 


5122905790 


006 






007 


ci 3cooonno 

J 1 J^UUUUUU 




000 


? f 4OO000OO 


5139637173 


009 


CI COOOOOOO 




ORDER 


COFFFICIFNT 


SIM OF SnilARF^ 


no 


ci t(^tL.(s\t\L. 






cn? 1 7 An? 1I1 

?U1 1 fOVJ It 


C9 t 399L9&C 


02 


4786911376 


l»8A4719fi1Q 


03 


4715197C9M 


RUtl 


TRANSFORMED Y 


PREDICTION 


001 


5111000000 


51 10066 5 82 


002 


5115000000 


5114543369 
5118504484 


003 


51 17000000 


004 


5120000000 


5120568866 


005 


5121*000000 


5122676974 
5129172463 


006 


5130000000 


007 


5135000000 


5135733204 


008 


51^0000000 


5140021700 


009 


5150000000 


5149912353 


OEFINEO X 


PREDICTION 




5150000000 


5118534484 




5210000000 


5129172463 




5215000000 


5140021700 







Fa^c j 



F RATIO x - 5188888889 
5170COOCCO 
5414324598 

F RATIO X - 5188888889 
516C000005 
5412317772 
4919324126 
RES ID. /ERROR 

5039623094 5110000000 
5033443777 5130000000 
511 56 603 6M 5150000000 
507389455M 5160000000 
5096961178 5170000000 
5072726522 5210000000 
503763203M 5213000000 
5031251165 5215000000 
501344837M 5220000000 

F RATIO X - 5188888889 
5150000001 
5411411195 
4917901880 
5055840582 
RESID. /ERROR 

4991710950 5110000000 
5037674665 5130000000 
51 12443 2 5M 5150000000 
50468 54 1 7M 51 60000000 
5110917863 5170000000 
5068781324 5210000000 
5057773 74M 5213000000 
491654702M 5215000000 
4957820963 5220000000 
ERROR ON Y 
5012090763 
5012031420 
5013114138 
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01010* 
01020 
01030 
01040XYZA 
01050 
01060XYZB 
01070 
01080 
01090 
01100XYZC 

01110 

01120 
01130 
01140 
01150 
01160 
01170 
01180A1 
01190M2 
012O0M3 
01210M4 
01220A5 
O1230A6 
01240A7 
01250A8 
01260A9 
01270 
01280A10 
01290 
01300A11 
01310 
01320A12 
O1330XYZ 
01340 
01350 



POLYNOMIAL REGRESSION 
DORG2178 
NOP ,6,10 
TF XYZB+6.XYZ 
AM *-1,10,10 
WATY 
RCTY 

TF XYZC+6.XYZ+5 
AM *-1,10,10 

WATY 
RCTY 
RCTY 

SM XYZA-1,1,10 
BNZ XYZA 
RCTY 
RNCD0 



B 

DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 



30, POLYNOMIAL REGRESSION PROGRAM? 
1+1, WRITTEN BY GENERAL FOODS RESEARCH CENTER? 



SWITCH 1 



SWITCH 2 
SWITCH 3 
SWITCH 4 



ON 
OFF 
ON 
OFF 
ON 
OFF 
ON 



DO RESIDUALS? 

NO RESIDUALS? 
DO PREDICTIONS? 

NO PREDICTIONS? 
PUNCH COEFFICIENTS? 

PRINT COEFFICIENTS? 
PUNCH RESIOUALS (SWITCH 1 
AND PREDICTIONS (SWITCH 2 ON)? 

OFF PRINT RESIDUALS (SWITCH 
A NO PREDICTIONS (SWITCH 2 ON)? 
40, PUT ONE CONTROL CARD AND DATA CARDS INTO 
15, READER HOPPER? 
26, PRESS READER START, START? 
DSA A 1 ,M2 ,M3 ,M4 ,A5,A6,A7,A8,A9,A 10 
DSA A11.A12 
DENO2176 



28, 
29, 
30, 

3 .'» 
34, 

35, 

Mt, 

31. 

hS, 

31, 



T 



ON) 
1 ON) 



o 



yj 



o 
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POLYNOMIAL REGRESSION PROGRAM 

WRITTEN BY GENERAL FOODS RESEARCH CENTER 

SWITCH f ON DO RESIDUALS 
OFF NO RESIDUALS 

SWITCH 2 ON DO PREDICTIONS 
OFF NO PREDICTIONS 

SWITCH 3 ON PUNCH COEFFICIENTS 
OFF PRINT COEFFICIENTS 

SWITCH if ON PUNCH RESIDUALS (SWITCH 1 ON) AND PREDICTIONS (SWITCH 2 ON) 
OFF PRINT RESIDUALS (SWITCH 1 ON) AND PREDICTIONS (SWITCH 2 ON) 

PUT ONE CONTROL CARD AND DATA CARDS INTO READER HOPPER 
PRESS READER START, START 



00903000000203 

1.11 

1.5 3 

1.7 5 

2.0 6 

2.4 7 
3.0 10 

3.5 13 
4.0 15 
5.0 20 
5 5 15 
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Y VERSUS X MINI'S 


8.888888 


B 


S 


2.688080 


SSE= 


a 1 


= 


.200950 


SS = 


B o 




2.603922 


SSE= 


e i 


: 


.208519 


SS = 


V. 2 


= 


.00011*7 


SS = 


■U'i 




TRANSFORMED Y 


PREDICTION 


001 




1.100000 


1.048122 


002 




1 .500000 


1 .461093 


003 




1 .700000 


1.075246 


00ft 




2.000000 


2.082765 


005 




2.1+00000 


2.290579 


006 




3.000000 


2.915792 


007 




3.500000 


3.543661 


OOO 




If. 000000 


3.963717 


009 




5.000000 


5.019021 


B 




2.6/4600 


SSE= 


B I 




.217603 


SS = 


B 2 




.000069 


SS = 


B 3 


SB 


.000151- 


SS = 


RUN 




TRANS FORMED Y 


PREDICTION 



OOJ 

002 
003 
004 
00 5 
006 
007 
003 
009 



.064623 
13.224265 



1. 100000 

1.500000 
1 .700000 

2.000000 
2.400000 
3.000000 
3.500000 
4.000000 
5 .000000 



OEFINED X 

5.000000 
10.000000 
15.000000 

nfADY FOR NF vr SET 



DF 
F 



1.006658 
1 .454336 
1.850448 
2.056006 
2.267697 
2.917246 
3.573320 
4.002170 
4.991235 

PREDICTION 

1.050440 
2.917246 
4.002170 



.064415 DF = 
13.224265 F = 
.000207 F . 



RES ID. /ERROR 

.396230 
.334437 
1.566036- 
.730945- 
.969611 
.727265 
.376320- 
.312511 
.134403- 



.057944 DF 

13.224265 F 

.000207 F 

.006471 F 



RES ID. /ERROR 

.091710 
.376746 

1.244325- 
.468541- 

1.091786 
.687813 
.577737- 
.016547- 
.057820 

ERROR ON Y 

.120907 
.120314 
.131141 



7.000000 
1432.459800 



6.000000 
1231.777200 
.019324 



5.000000 
1141.119500 
.017901 
.558405 



o 
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020!0 
0202C 

02030 DS 
020U00UTPUT0S 
02050 INPUT DS 
020601 DS 
020700UT1 DS 
O2OB0OUT2 
020SOOUT3 
02100 
02110T1 
02t20T2 
02130T3 
02140T4 
02150T5 
02160T6 
02170T7 
02100Y 
021 90X 
02200N 
022100 

02220REVYX DS 
02230TP.Y DS 
02240TRX DS 
022SOSTARTRDS 
02260STARTPDS 
02270CR0SS DSB 
022806 DS 
02290T0TAL DS 
023000FE DS 
02310SUMS DSB 
02320VECTORDS 
02330SSB DS 
02340TOTSS DS 



DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 
DS 



02350SSE 
02360MM 
02370AVGY 
02380MSE 
02390AVGX 
24000 VAR 



02410FLZER00C 



0242081 
0243082 
02440B3 
021*5064 
02460644 
0247065 
02480655 
0249066 
0250068 



SYMBOLS AND CONSTANTS 
D0RG2178 
2 
1 
1 
2 

15 
15 
15 
111 
10 
10 
10 
10 
10 
10 
10 
10 
10 
3 
2 
1 
2 
2 
2 
2 

10,10 
I, CROSS 
10 

1 .TOTAL 
10.19 
1,SUMS 
1,SUMS+110 
10 

1.SUMS+100 
10,55 
10 

1 ,AVGY 
10 

2 

10,0 

10, Y VERSUS 
9,X MINUS 
7,6 1 - 
8, SS - 
8, SSE- 

7. £ - * 
OF - 



DS 
DS8 
DS 
DS 
DS 
DS 



DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 
DAC 

02510HEA01 DAC 49, ORDER 
02520HEA01ADAC 7, 1 1 
02530HEA018DAC 5, X ♦ 
02540HEAD1CDAC 20, 1111111111 
02550HEA02 DAC 50, RUN TRANSFORMED Y 
02560 DAC 31, 

0257CHCAD3 DAC 47, DEFINED X 
025PO DAC 34, 

025S0R0UHD DC 10,1 



7, OF - 

4 1110 

19, READY FOR NEXT SET0 

COEFFICIENT SUM OF SQUARES 



F RATIO 



PREDICTION 


PREDICTION 



RESID./ERROR0 
ERROR ON Y0 
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03C 10* TRANSFORMATION SUBROUTINE 

03020 DAS 5 

03030TRANS B 

03040 D0RG*-3 

03050LOG FL0GTRANS-1 , TRANS-1 

03060 B AS I S 

03070 D0RG*-3 

03080EXP FEX TRANS-1 ,TRANS-1 

03090 B AS IS 

03S 00 D0RG*-3 

03110RECIP FD FLONE, TRANS-1 

03120 TF TRANS-1, 99 

03130ASIS BB 

03140 D0RG*-9 

03150ADDR DS 5.TRANS+6 

04010* SET UP TRANSFORMATION BRANCH. AODR CONTAINS LABEL TO BRANCH TO. 

04020 DAS 1 

04030SET TFM AD DR, AS IS 

04040 CM SET-1,1,10 

04050 BN SET3 

04060 BH SET1 

04070 TFM ADDR.LOG 

04080 WATYA2 

04090 B SET 3 

04100 D0RG*-3 

04110SET1 CM SET-1,2,10 

04120 BH SET2 

04130 TFM ADDR.EXP 

04140 WATYA3 

04150 B SET3 

04160 D0RG*-3 

04170SET2 TFM ADDR,REC I P 

04180 WATYA4 

04190SET3 BB 

04200 D0RG*-9 

0421 OA 2 DAC 14, LOGARITHM OF 

04220A3 DAC 10, E TO THE 

04230A4 DAC 15,REC IPROCAL OF 

05010* MATRIX ADDRESS COMPUTATION. 99 IS DESTROYED. 



05020V1 


DAS 


1 


05030V2 


DAS 


1 


05040MAC 


TFM 


AODR, MM 


05050 


TF 


V3.V1 


05060 


C 


V2.V1 
MAC1 


05070 


BNN 


05080 


TF 


V3,V2 


05090 


TF 


V2.V1 


05100MAC1 


A 


ADDR-1 ,V2 
V2,19,10 


05110 


TFM 


05120 


SM 


V2 


05130 


MM 


V3,5,10 


05140 


SF 


98 


05150 


TF 


*+23,99 


05160 


MM 


V2 


05170 


A 


ADD R, 99 


05 180 


BB 


05190 


D0RG*-9 



Pegs h 



FA 
AM 

TF 



06100W0RK2 FM 
06110 
06120 
06130 
06140W0RK3 BT 
06150 TF 
06166 FM 
06170W0RK4 FM 



05200V3 OS 1, MAC 1+35 

06010* WORK SUBROUTINE USED FOR RESIDUALS AND PREDICTIONS 
06020WORK FS X,AVGX 
06030 TF T1.FL0NE 

06040 TFM WORK2+23.B 

06050 TFM J, 0,10 

06060 TF 0UT2.FLZER0 

06070 TF T2.FL0NE 

06080W0RK1 TF Vl.J 
06090 TF T3.T1 

B.fl 
OUT 2, 99 
WORK2+23 ,10,10 

K,J 
MAC,K 

W0RK4+35.ADDR 
T1.T3 
99.MM 

06180 FA T2,99 

06190 C J.K 

06200 BE WORKS 

06210 FA T2,99 

06220W0RK5 FM T3,X 
06230 TF T3.99 
06240 AM K.I. 10 

06250 C K.DVAR 

06260 BNH W0RK3 

06270 FM T1.X 

06280 TF T1.99 
06290 AM J, 1,10 

06300 - C J.OVAR 
06310 BNH WORM 

06320 FM T2.MSE 

06330 TF T2.99 

06340 FSQR0UT3.T2 
06350W0RK7 B ODD 
06360 D0RG*-4 
06370 DNB 50 

06380 DNB 30 

06390 DC 1 ,9 

06400BLANK OS 1.WORK7+7 

07010* BEGINNING OF MAIN PROGRAM 

07020START RCTY 

07030 RNCDN-2 

07040 SF N-2 

07050 SF D-1 

07060 SF TRY-1 

07070 SF TRX-1 

07080 SF STARTR-1 

07090 SF STARTP-1 

08010* SET BRANCHES FOR Y AND X TRANSFORMATIONS 

08020 BT SET, TRY 

08030 TF ABB+1 1 ,ADDR 

08040 WATY B1 

08050 BT SET,TRX 

08060 TF ACC+11,ADDR 

08070 WATY B2 
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09010* CLEAR AND SET UP ACCUMULATORS 

09020 TFM CLEAR+6, CROSS 

09030CLEAR TF ,FLZERO 

09040 AM *-6,10,10 

09050 CM CLEAR+6.CROSS+330 

09060 BN CLEAR 

10010* DATA INPUT LOOP 

10020 TF CI.N 

10030 TFM AC+6.DATA 

10040 TFM AD+6.DATA+10 

10050AB TFM CCCC2+6,Y,,SET STORE 

10060 BTM CCCC,2,10,TW0 INPUTS 

10070 BD DOREV.REVYX 

10080 B ABB 

10090 D0RG*-3 

IfllOODOREV TF T5,X 

10110 TF X,Y 

10120 TF Y,T5 

10130ABB TFM ADDR 

10140 BT TRANS, Y 

10150AC TF .TRANS-1,, STORE TRANSFORMED Y 

10160 AM *-6, 20,10 

10170 FA MM, TRANS-! 

10180ACC TFM ADDR 

10190 BT TRANS, X 

10200AD TF ,TRANS-1,, STORE TRANSFORMED X 

10210 AM *-6,20,!0 

10220 FA MM+10.TRANS-1 

10230 FA TOTAL, FLONE 

10240 SM CI, 1,10 

10250 BNZ AB 

11010* GET AVERAGES. COMPLETE STATEMENT ON COOING OF X 

11020 FD MM+IO, TOTAL 

11030 TF AVGX, 99,, STORE AVERAGE X 

11040 TF ROUND-8,91 

11050 FA AVGX, ROUND 

11 060 AD 2 BT WWWW.AVGX, .PRINT AVERAGE X 

11070 FD MM, TOTAL 

11080 TF AVGY,99, STORE AVERAGE Y 

11090 RCTY 

12010* READ X LOWER, DELTA X, AND X UPPER FOR PREDICTIONS 

12020 TDM EEI+6,0 

12030 BNC2AD4 

12040 TDM EE 1+6,1 

12050AD3 TFM CCCC2+6,T5,,SET STORAGE LOCATION 

12060 BTM CCCC,3.10,READ X LOWER, DELTA, UPPER 

12070 BD A04.EE1+6 

12080 TDM EEI+6,1 

12090 B EE1 

12100 D0RG*-3 

13010* SUM POWERS OF X. SUM CROSS-PROOUCTS OF Y AND POWERS OF X 

13020AD4 TF C1,N 

13030 TFM AF2+11.SUMS+10 

13040 A AF2+10,D 

13050 A AF2+10.D 



o 



13060 


TFM 


AH2+1 1 ,CR0SS+10 


13070 


A 


AH2+10.D 


13000 


TFM 


AE+11 ,DATA+10 


13090 


TFM 


AG+1 1 ,DATA 


131 00AE 


TF 


X,,,PICK UP X 
*-1,20,10 


13110 


AM 


13120 


FS 


X.AVGX 
AF+23,SUMS 


13130 


TFM 


13140 


TF 


99, FLONE 


13150AF 


FA 


SUMS, 99 


13160 


AM 


AF+23. 10. 10 


131 70 
13J80AF2 


FM 




CM 


AF+23 


13190 


BN 


AF 


13200 


TFM 


AH+23, CROSS 


13210AG 


TF 


Y,,,PICK UP Y 


13220 


AM 


*-1,20,10 


13230 


TF 


99, Y 


13240AH 


FA 


CROSS, 99 


13250 


AM 


AH+23, 10, 10 


13260 


FM 


99.X 


13270AH2 


CM 


AH+23 


13200 


BN 


AH 


13290 


FS 


Y.AVGY 


13300 


FM 


Y,Y 


13310 


FA 


T0TSS,99 


13320 
13330 


SM 


CI, 1,10 


BNZ 


AE 



14010* SET 

14020 

14030 

14040 

14O50AJ 

14060 

14070AK 

14080 

14090 

14100 

14110 

14120 

14130 

14140 

15010* SET 

15020 

15030 

15040 

15050 

15060 

15070 

15080 

15090 

15100 

151 fO 

15120 

15130 

15140 

15150BBB1 

15160 



UP X PRIME X MATRIX 

TFM AJ+11 ,SUMS, , INITIALIZE PICKUP 

TFM AK+6,MM,, INITIALIZE STORE 

TFM CI, 10, 10, SET BASIC COUNTER 

TFM AK+11,,,SET PICKUP FOR DIAGONAL 

TF C2,C1,,SET SECONDARY COUNTER 



AM AK+6,10,10 
AM AK+11,10,10 
SM C2,l,10 
BNZ AK 

AM AJ+11, 20, 10 
SM CI, 1,10 
BNZ AJ 

UP LOOP FOR INCREASING ORDER OF POLYNOMIAL 

TFM DVAR, 1,10, LIMIT IS INPUT D 

FD CROSS, MM 

TF B, 99,. SET FIRST B 

TFM BD+11.CR0SS+10 

TF SSE,TOTSS 

TFM BN+6,SSB 

FS TOTAL, FLONE 

FD FLONE, MM 

TF MM, 99,, SET 1 BY 1 INVERSE 

TD HEAD 1A+6, TRY 

TD HEAD1A+12.TRX 

TFM BBS 1+6, HEAD 1C+2 

TFM BBB1+11,AVGX-9 

TD 

AM BBB1+6.2.10 



15170 AM BBB 1+1 1,1,10 

15180 CM BBB1+11.AVGX+1 

15190 BN BBB1 

15200 TFM HEADIB+8,20,10 

15210 BNF *+24,AVGX 

15220 TFM HEADIB+8,10,10 

15230BBBBBBBNC3BB-24 

15240 TR OUTPUT, BLANK 

15250 WACDHEAD1 

15260 B BB 

15270 D0RG*-3 

16010* GET AUXILIARY VECTOR, START NEW DIAGONAL ANO NEW B 



16020 


RCTY 


1 £nan 


RCTY 


lOUtUDD 


TFM 


1,0,10 


1 o05O 


TF 


VI, DVAR 


1 andn 


BT 


MAC, DVAR 


1 OU/U 


TF 


BC+11,ADDR,,DIAG 


1 Anftn 


TF 


BL+6.ADDR, ,D IAG 


1 ouyuDv. 


TF 


T1,MM,,PICK UP DIAG 


1 AinnRn 

1 O IWBV 


TF 


T2,B+10,,PICK UP XY CROSS-PRODUCT 


i A 1 in 

1 I IU 


TFM BK+23,B,,NEW B 


1 Ai on 


TFM BH+6, VECTOR 




TF 


T3,FLZER0, , VECTOR 


i£iLn 

ID IHU 


TFM 


J, 0,10 


1 A 1 cnR r 


TF 


VI, 1 


1 £.1 en 


BT 


MAC, J 

BG+23 ,ADDR 


lol /U 


TF 


\ D 1 ou 


TF 


VI, J 


1 A ion 

IOI J\J 


BT 


MAC, DVAR 
BG+35.ADDR 


1 OZ.W 


TF 


1 A 9 inR£ 
1 Of. Iviuu 


FM 


MM, MM 
T3.99 


1 (.ton 


FA 




AM 


J, 1,10 


16240 


C 


J, DVAR 


16250 


BN 


BF 


16260BH 


TF 


VECTOR, T3,, STORE AUX VECTOR 


16270 


AM 


*-6,10,10 


16200 


TF 


VI, 1 


16290 


BT 


MAC, DVAR 
BJ+35,ADDR 


16300 


TF 


16310 


TF 


BK+35,ADDR 


16320BJ 


FM 


T3.MM 


16330 


FS 


T1,99 


16340BK 


FM 


B,MM 

BK+23, 10,10 


16350 


AM 


16360 


FS 


T2.99 


16370 


AM 


1,1.10 


16300 


C 


1 ,DVAR 


16390 


BN 


BE 


16400 


FD 


FLONE, T1 

MM, 99,, STORE NEW DIAGONAL 


16410BL 


TF 


16420 


TF 


T3, 99,, SAVE NEW DIAGONAL 


16430 


TF 


T4,99 


16440 


SF 


T4 


16450 


TF 


BM+6.BD+11 


16460 


AM 


BD+11,10,10 
T2.T3 


16470 


FM 


16400BM 


TF 


B+10, 99,, STORE NEW B 



161*90 FM 99.T2 

16500BN TF SSB, 99,. STORE SUM OF SQUARES 

16510 AM *-6,10,lO 

17010* FORM ERROR D.F. AND SUM OF SQUARES 

17020 FS SSE,99 

17030 FS DFE,FLONE 

17040 FO SSE,DFE 

17050 TF MSE, 99,, STORE ERROR VARIANCE 

18010* ENLARGE INVERSE, CORRECT PRIOR INVERSE, OUTPUT 

18020 TFM CJ+23,B 

18030 TFM CE+23, VECTOR 

18040 TFM CLL+11,SSE 

18050 TFM 1,0,10 



19010* CORRECT 
19020CC TF 
19030 
19C40 
19050CD 
19060 
19070CE 
19080 
19090 
19100CF 
19110 
19120 
19130 



TF 
TF 
BT 
TF 
FM 
AM 
FM 
FA 
AM 
C 

BN 



PRIOR INVERSE 
J, I 

CE+35.CE+23 
VI, I 
MAC, J 

CF+23,ADDR 

VECTOR, VECTOR 

CE+35,10,10 

99.T3 

MM, 99 

J, 1,10 

J,OVAR 

CO 



20010* ENLARGE INVERSE 

20020 BT MAC,OVAR 

70030 TF CH+6,AD0R 

20040 TF CG+23, CE+23 

20050 AM CE+23, 10. 10 

20060CG FM VECTOR, T4 

20070CH TF MM, 99,, STORE NEW COLUMN OF INVERSE 

21010* CORRECT PRIOR B VALUES 

21020 FM 99, T2 

21030CJ FA B,99 

21040CK TF CL+11,CJ+23 

21050CL TF 0UT1,B,,PUT B INTO OUTPUT AREA 

21060 AM CJ+23,10,10 

21070CLL TF 0UT2,SSE,,PUT SUM OF SQUARES INTO OUTPUT AREA 

21080 AM *-l,10.10 

21090 FD 0UT2,MSE 

21100 TF 0UT3;99,,PUT F RATIO INTO OUTPUT AREA 

21110 BNC3CM 

21120 CF 0UT1-9 

21130 CF 0UT2-9 

21140 CF OUT3-9 

21150 CF 1-1 

21160 BNF CLL2.00T! 

21170 TDM 0UT1,4,11 

21180CLL2 WNCDOUTPUT 

21190 SF 1-1 

21200 B CN 

21210 D0RG*-3 

21220CM RCTY 



e 



21230 


TD B3+4.I 


21240 


WATYB3 


21250 


BT WWWW.0UT1 


21260 


CM 1,0,10 


21270 


BNE CMt-12 


21280 


WATYB44 


21290 


B CM1 


21300 


D0RG*-3 


21310 


WATYB4 


21320CM1 


BT WWWW,0UT2 


21330 


CM 1,0,10 


21340 


BNE CM2-12 


21350 


WATYB55 


21360 


B CM2 


21370 


D0RG*-3 


21380 


WATYB5 


21390CM2 


BT WV/V/W,0UT3 


21400CN 


AM 1,1,10 


21410 


C 1 ,OVAR 


21420 


BN CC 


21430 


TF V1,DVAR 


21440 


BE CK 


21450 


RCTY 


21460 


C DVAR.STARTR 


21470 


BN EE-24 


21480 


BNC1 EE-24 



22010* RESIDUALS PROGRAM 

22020 TFM W0RK7+6.D0D 

22030 TFM DD+23 ,DATA 

22040 TFM DD+35,DATA+10 

22050 TFM 1,1,9 

22060 BNC4DD1 

22070 WACDHEAD2 

22080 B DD 

22090 D0RG*-3 

22100DD1 RCTY 

22110 RCTY 

22120 WATYHEAD2 

22130 RCTY 

22140DD BNC1 EE-24 

22150 TF 0UT1,,,PICK UP Y 

22160 TF X 

22170 TF OUT3+12,X 

22180 AM DD+23, 20 ,10 

22190 AM DD+35,20,10 

22200 B WORK 

22210 D0RG*-3 

222200DD TF 99,0UT1 

22230 FS 99,0UT2 

22240 FD 99.0UT3 

22250 TF 0UT3, 99,, STORE STANDARD I ZED RESIDUAL 

22260 BNC4DE 

22270 CF OUT 1-9 

22280 CF OUT2-9 

22290 CF OUT3-9 

22300 CF OUT3+3 

22310 CF 1-2 

22320 BNF DDD1.0UT1 

22330 TDM OUT! ,4,11 



o 



2234C0DD1 

22350 

223600002 

22370 

2238ODO03 

22390 

221*00 

221*10 

22420DE 

22430 

22440 

22450 

22460 

22470 

22480 

22490 

225O0DF 

22510 

22520 



BNF DD02,OUT2 
TOM OUT 2, 4, 11 
BNF OD03.0UT3 
TOM OUT3.4.1I 
WNCDOUTPUT 
SF 1-2 
B OF 
D0RG*-3 
TO B6,l-2 
TO B 6+2, 1-1 
TO 86+4,1 
RCTY 
WATYB6 

BT WWWW.0UT1 
WWWW.0UT2 
WWWW.0UT3 
».1.10 
I.N 



BT 
BT 
AM 

C 



BMH 00 



23010* PREDICTIONS PROGRAM 
23020 
23030 
23040EE 
23050 
23060 
23070 
23080 
23090 
23100 
23110 
23120N0G0 
23I30EE2 



C DVAR.STARTP 
BN FF 
BNC2FF 

BO EE2.EE1+6 
RCTY 
RCTY 

WATYNOGO 
RCTY 
B AD 3 
D0RG*-3 

DAC 3 2, READ PREDICTION DEFINITIONS NOW? 
BNC4EE1 



23140 TR OUTPUT, BUNK 

23150 WAC0HEA03 

23160 B EF-24 

23170 0ORG*-3 
23180EE1 RCTY 

23190 RCTY 

23200 MATYHEA03 

23210 RCTY 

23220 TF T4.T5,, INITIALIZE X LOWER 

23230 TFM WORK7+6.EG 

23240EF BNC2FF 

23250 TF AOOR.ACC+11 

23260 BT TRANS, T4,, TRANSFORM X 

23270 TF X.TRANS-1 

23280 B WORK 

23290 D0RG*-3 

23300EG BNC4EF1 

23310 TF 0UT1.T4 

23320 CF OUT 1-9 

23330 CF OUT 2-9 

23340 CF 0UT3-9 

23350 BNF *+24,0UT1 

23360 TOM OUT 1,4, 11 

23370 BNF *+24,OUT2 

23380 TDM OUT 2, 4, 11 

23390 WNCDOUTPUT 

23400 B EG 2 



© 




23410 D0RG*-3 
23420EF1 RCTY 

23430 BT WWWW,T4, .PRINT X 

23440 BT WWWW.0UT2, .OUTPUT PREDICTION 

23450 BT WWWW,0UT3,,0UTPUT STD. ERROR 

23460EG2 FA T4.T6 

23470 C T4.T7 

23480 BNH EF 

23490FF AM DVAR,1,|0 * 

23500 C DVAR.D 

23510 BNH BBBBBB 

23520 RCTY 

23530 WATYB8 

23540FG H 

23550 B START 

23560 D0RG*-3 

23570J DS 2,START+35 

23580K DS 2,START+47 

23590C1 DS 5.START+59 

23600C2 DS 5.START+71 

23610FL0NE DC 10,51 10000000, FG+1 1 

24010* SUBROUTINE FOR FLOAT TO FIXED TO PRINT 484 DIGITS 

2402CWWWWA DAC 8,00000000 

24030 DAC 8,000000 <3> 

24040WWWWB DAS 1 

24050WWWWD DS 16 

24060WWWWE DS 12 

24070 DS 4 

24080WWWWC DS 2 

24090 DS 8 

24100WWWW TR WWWV/B-1 .WWWWA-1 

24110 CM WWWWC.59.10 

24120 BN WWWW1 

24130 S WVAA/E,WWWWE 

24140 WATYWWWWD+2 

24150 TF WWWWE+1,WWWW-1 

24160 WNTYWWWWD+4 

24170 BB 

24180 D0RG*-9 

24190WWWW1 TFM WWWW3+11,WWWWC+1 

24200 TFM WWWW3+6,WWWWD 

24210 SM WWWWC,50,10 

24220 S WWWW3+6.WWWWC 

24230 S WWWW3+6.WWWWC 

24240 TD WWWW8+7,WWWW-1 

24250 CF WWV/W-1 

24260WWWW2 CM WWWW3+6.WWWWE 

24270 BNN WWWW4 

24280 CM WWWW3+1 1 ,WWWW 

24290 BE WWWW4 

24300WWWW3 TD 

24310 AM WWWW3+11,1,10 

24320 AM WWWW3+6,2,10 

24330 B WWWW2 

24340 D0RG*-3 
24350WWV/W4 TF WWWWE+4,WWWWE+2 
24360 TFM WWV/WD.3,10 

24370 BNF WWWW5.WWWW8+7 

24300 TFM VWWWE+2,20.10 



Yjw 11 



2439CWWWW5 

24400WWWW6 

24410 

24420WWWW7 

24430 

24440 

24450 

24460WWWW8 

24470 

24480 



TFM WWWW6+H.WWWWB 
BD WWWW8 

TF WWWW7+6.WWWW6+11 

TFM ,0,10 

AM WWWW6+11,2,10 

CM WWWW6+1 l.WWWWD 

BN WWWW6 

WATYWV/WWB 

BB 

D0RG*-9 



25010* CONVERSION OF VARIABLE LENGTH FIXED POINT DATA FROM CARDS 



25020 


DAS 


1 


25030CCCC 


TFM CCC 1+6, INPUT 


25040CCC1 


TFM 


.0,10 


25C 50 
25060 


AM 


CCCI+6,2,10 
CCC1+6.INPUT+144 


CM 


25070 


BNZ 


CCCf 


25080CCC2 


RACD INPUT 


25090 


TR 


INPUT+143.CCCCCC+1 


25100CCC3 


TR 


INPUT-3, INPUT-1 


25110 


BNR CCC4, INPUT -2 


25120 


B 


CCC2 


25130 


D0R6*-3 


25140CCC4 


CM 


INPUT-2,0,10 


25150 


BE 


CCC3 


25160 


TR 


CCCCH-l.CCCCCC-11 


25170 


B 


CCC6 


25180 


D0RG*-3 


25190CCC5 


TR 


INPUT-3, INPUT-1 


25200CCC6 


CM 


INPUT-2,70,10 


25210 


BN 


CCCC3 


25220 


TR 


CCCD3.CCCD3+1 


25230 


TD 


CCCDDD, 1 NPUT-2 


25240 


BD 


CCC7.CCCD2 


25250 


AM 


CCCCH',1,10 

CCC5.CCCD3 


25260CCC7 


BNF 


25270 


CM 


INPUT, 45, 10 


25280 


BNE 


CCC8 


25290 


TD 


CCC777+1.INPUT+1 


25300 


TD 


INPUT+5, INPUT+4 


25310 


SF 


1 NPUT+5 


25320CCC777A 


CCCCH.INPUT+6 


25330 


TR 


INPUT-3, INPUT+5 


25340CCC8 


C 


CCCDDD, FLZERO 


25350 


BE 


CCCC2 


25360CCC9 


TF 


CCCD3+1 ,CCCCH 


25370 


BD 


CCCC 1 ,CCCD3+2 


25380 


SM 


CCCCH,1,10 
CCCD3.CCCD3+1 


25390 


TR 


25400 


TDM 


CCCDDD.O 


25410 


B 


CCC9 


25420 


D0RG*-3 


25430CCCC1 


BD 


CCCC2,CCCD1 


25440 


SF 


CCCDDD 


25450CCCC2 


TF 


,CCCDDD 
CCCC2+6,10,10 


25460 


AM 


25470 


SM 


CCCC-1,1,10 


25480 


BNZ 


CCC3 


25490 


BB 





25500 D0RG*-9 
25510CCCC3 CM INPUT -2, 20, 10 
25520 BNE CCCC4 

25530 TDM CCCD1,0 

25540 B CCC 5 

25550 D0RG*-3 
25560CCCC4 CM INPUT -2, 3, 10 
25570 BNE CCCC5 

25580 TDM CCCD2, 1 

25590 B CCC5 

25600 D0RG*-3 
25610CCCC5 CM INPUT-2,0,10 
25620 BNE CCCC8 

25630CCCC7 TR CCCD3 ,CCCD3+1 
25640 TDM CCCDDD, 

25650CCCC6 BNF CCCC7,CCCD3 
25660 B CCC8 

25670 D0RG*-3 
25680CCCC8 CM I NPUT-2, 10, 10 
25690 BE CCC5 

25700 WATY I NPUT-2 

25710 H 
25720 D0RG*-9 
25730CCCCCCDC 12,501000000000 
25740 DC 2,0 

25750 DAC 1,(5) 

25760CCCCH DS 2 
25770CCCD1 OS f 
25780CCCD2 DS 1 
25790CCCD3 DS 1 
25000CCCDDDDS 9 
25810CCCD4 DS 2 

26010* 

26020DATA DSB 10,500 
26030 DEND2178 
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